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ABSTRACT 


A model forecasting the wave height and the longshore current distribution inside 
the surf zone based on the formulations by Thornton and Guza (1983), (1986), are ap- 
plied to an extensive set of both laboratory and field data for the purposes of testing and 
modification. The models are tested on planar beaches as well as barred beaches for a 
variety of wave conditions. 

The wave transformation model is based on solving the energy flux equation using 
a bore dissipation mechanism and describing the random wave heights with the Rayleigh 
distribution. The two model parameters B and y, where B describes the amount of foam 
of a breaking wave and jy is the proportionality constant which relates the rms wave 
height to the water depth, are combined into a single parameter BG. The combined 
parameter BG is shown to be a function of deep water surf similarity parameter. Ap- 
plied to the present data sets, the rms error of the measured wave height and the model 
predicted wave height was usually less than 9°%> and ranged from 1.5% to 15.7% with 
amean of 5.3°o anda standard deviation of 3.1% for the whole 74 data sets. The 
Wave transformation model is highly robust in describing the wave height distribution in 
the surf zone. 

The longshore current model is based on solving the steady state, alongshore mo- 
mentum balance for straight and parallel contours using the radiation stress concept. 
The model requires specifying the bed shear stress cocfficient (c,) and turbulent mixing 
coefficient (N). Applied to the present data sets the rms error between the measured and 
modeled longshore current values ranged from 4.5% to 55.5% witha mean of 24.6%. 
The turbulent mixing is not required for planar beaches, but it is required for barred 
beaches to describe the longshore currents inside the surf zone. The mean value of the 
bed shear stress coefficient for the present longshore current data is 0.006 and the mean 
lateral mixing coefficient is 0.006. 
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I. INTRODUCTION 
A. BACKGROUND 


The wave height and longshore current distribution inside the surf zone is a primary 
concern of Coastal Engineering and Naval Operations. These distributions are impor- 
tant for amphibious operations both at the decision making stage and during the beach 
assault in order to reduce the percentage of casualties. The type of landing craft to be 
chosen depends on the wave and current conditions in the surf zone. 

As waves move from deep water to shallow water, they are transformed due to the 
effects of shoaling. refraction and dissipation. Shoaling results in decrease in wave 
length and an increase in wave height. Wave refraction results in change in wave direc- 
tion. Shoaling and refraction are reasonably well defined by theory. Dissipation is due 
to bottom friction and wave breaking. Wave breaking is the dominant dissipative 
mechanism inside the surf zone. At a certain depth the waves become unstable and 
break in either spilling. plunging or surging form, depending on the deep water wave 
steepness and the beach slope. After breaking the wave height decreases as the energy 
is converted into turbulence and finally lost to heat. If the waves approach the beach 
at an angle then thev cause a longshore current parallel to the shoreline within the surf 
zone because of the change in the alongshore component of momentum. 

The motion of water inside the surf zone is extremely complex because of the dy- 
nanuc upper and lower boundaries and generally three dimensional nature of the flow. 
In order to overcome this complexity in the surf zone, a number of assumptions have 
to be made. As a consequence. the models describing wave transformation and 
longshore currents include empirical parameters that must be specified in the formu- 
lations. The specification of these empirical parameters requires testing the model 
against comprehensive experimental data. 

As of today, there is no exact theoretical formulation for the wave height distrib- 
ution inside the surf zone where the waves break. Therefore empirical formulations are 
used to specifv the wave height in the surf zone. It is generally assumed that the 


breaking wave height in the surf zone is a linear function of the depth in the form of: 


I] = yh (1) 


in which H is the wave height, h is the water depth and y is the proportionality constant. 
This relation holds fairly well for spilling breakers. It has been shown by several inves- 
tigators such as Horikawa and Kuo (1967), Nakamura, Shirashi and Sasaki (1967), 
Divorky, Le mehauté and Lin (1970) that this relation does not hold for plunging 
breakers. Although this relation does not hold for all breaker types it has been retained 
in most models built to predict the wave height in the surf zone (Smith and Kraus, 1988). 
Wu_ et al. (1985) also utilize this relation in their two dimensional model that includes 
the bottom friction and the nonlinear convective accelerations. The major deficiency of 
present day models is the inability to correctly describe wave breaking processes. 


B. OBJECTIVES 

The objectives of this thesis are to test and modify the wave transformation and 
longshore current models being considered for application in the Tactical Environmental 
Support System (TESS). These models are part of the the surf forecasting software 
which was developed under contract by Marshall D. Earle, MEC Systems Corporation 
to provide the Navy with an improved Surf Forecasting capability. The basis of this 
model are the wave transformation and longshore current models developed by 
Thornton and Guza (1983), (1986). In order to make the models more operational and 
increase the range of applicability. an expanded data base of 74 different data sets in- 
cluding field and lab data are utilized. The model parameters are calculated by mini- 
mizing the error in the least square sense between the model and the data, an attempt 


is made to predict the model parameters. 


N 


Il. NUMERICAL MODEL 


The surf forecasting model predicts the wave height and longshore current distrib- 
ution across the surf zone. The feature of this model compared to most previous studies 
is that the waves are treated as random rather than monochromatic. The coordinate 
system orientation and sign convention used in the model is given in Figure 1. The 
shoreline has been aligned with the v-coordinate and the x-coordinate is positive off- 
shore. A wave orthogonal approaching the beach from southerly direction has a positive 
angle value (+@) and causes positive longshore current, ie, to the right, where as a 


northerly wave orthogonal is just the opposite. 


A. DESCRIPTION OF WAVE TRANSFORMATION MODEL 

The wave height transformation for straight and parallel contours is calculated using 
the energy flux balance in which the onshore spatial variation of the energy flux is bal- 
anced by the ensemble averaged dissipation due to wave breaking, < e, > , and the 


frictional dissipation. < e, > . The energy flux balance is given by: 


ere.) 


=<6>+<e-> p. 
Px b f (2) 


where E is the wave energy, C,, 1s the x-component of the group velocity. 

Based on analvsis of field data acquired at Torrey Pines Beach, CA, it was found 
that the Rayleigh distribution was in good agreement with the data both outside and 
inside the surf zone. Therefore the wave height distribution inside the surf zone is de- 


scribed by the Rayleigh probability density function given by: 





pt) =o exp( ~( a [3] 


rms 


Energy and the x-component of the group velocity are described using linear theory. 


r= + pg H p(H)dil = = ps Hems 


ry 





Figure 1. Orientation of the Coordinate System and Sign Convention: The y-axis 


has been ahgned with the alongshore and the x-distance is positive off- 
shore. A wave orthogonal coming from the south has a positive angle 
value and causes a longshore current to the nght as the landing craft 


approaches from offshore to the beach. 


Qkh 
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On = > es 
where k is the wave number associated with peak frequency /, h is the water depth and 
0 is the mean wave direction corresponding to the peak of the wave spectrum. 
Although there are several dissipation mechanisms (frictional, percolation, sediment 


transport, etc. ) once the waves start to break turbulent dissipation becomes the domi- 


nant dissipative mechanism. Therefore the model only includes dissipation due to wave 
breaking. The frictional dissipation is not included in the model because it was found 
that the frictional dissipation was less than 3°% of the dissipation due to wave breaking. 
Even if this is not always true, the effect is compensated for by the B parameter in the 
bore dissipation mechanism. It is more reasonable to include the frictional effects in B 
rather than introducing new parameters. Frictional effects may become significant in the 
swash regions but this is not the concern of this model. 

Unlike monochoromatic waves in which there is a single breaker line, the random 
treatment of the waves describes the larger waves to break farther offshore and the 
smaller waves to break closer to the shore line, which is more realistic. Across the surf 
zone at each point there might be broken and unbroken waves at the same time. The 
model calculates the breaking wave dissipation utilizing the periodic bore formulation 
as suggested bv LeMehauté(1962). The bore dissipation mechanism as applied only to 


the breaking waves is formulated as follows : 


4 (y-hy 1 (Bye 
Ebore = et Ps at ees a Ps 


[6] 
where p is the water density, g is the gravitational acceleration, A,1s (h +) , A, is 
(h-» ) , 9 1s the surface elevation, h is the mean water depth and Q is the volume 
discharge across the bore. Measurements done on a gently sloping (1:40) beach show 
that the classical periodic bore formulation underestimates the dissipation by 30% to 
S0°o (Stive 1983). The under-estimation is compensated for by adjusting the value of 
the parameter B which describes the amount of foam on the crest of a wave (see 
Figure 2). There are several formulations suggested for the volume discharge across the 


bore. The description by Hwang ez al. (1970 ) is used in the model, 


a (7) 


where C is the wave speed, L is the wave length and h is the mean water depth. 
In order to find the percentage of waves that are breaking, the probability density 


function is simply multiplied by a weighting funtion. 


PID) = WUD pH) [8] 


tn 





Figure 2. Sketch of the Periodic Bore Formulation: [1 is the wave height, h is the 


mean Water level, /, its the water depth below wave trough, /, is the 
surface elevation measured from the bottom. B is the amount of foam 


over the crest of a breaking wave. 


Note that 7,(/7), 1s not a pdf but a percentage. In order to distinguish the percentage 
of the waves that are breaking, the following empirical weighing function which fits the 


observations is used (Thornton and Guza, 1983) 


i] 
WUD = Cat — exp -( 4) a [9] 
It was found by the same authors that a value of y=.40 and n=2 gave the best fit to 
the data. 
The average rate of energy dissipation is found by adding up the dissipation for each 
breaking wave calculated using the formula 


(BIN 


fr 
i aa [10] 


[bore = 


Then this sum ts divided by the total number of waves both broken and unbroken. In 


Other words the average rate of energy dissipation is calculated by multiplying the dissi- 


pation for a single broken wave height H by the probability of wave breaking at each 
height, as given bv P,(H). For the ensemble 


S B: i > 3 
o> =—7 ps | HP (H)dH C11] 
0 


Substituting (10) the bore dissipation is given by : 
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B. DESCRIPTION OF THE LONGSHORE CURRENT MODEL 
Most of the longshore current models for straight and parallel contour beaches are 
based on an alongshore momentum balance. Assuming steady state, the alongshore 


momentum balance is given: 








cs, eS 
ee ee C13] 


in which the first term on the left hand side is the wave induced momentum flux, or ra- 
diation stress gradient and the second term is the lateral transfer of turbulent momen- 
tum. On right hand side of the equation, are the bottom shear and surface wind stresses 
in the alongshore direction. 

The wave induced radiation stress due to gravity waves, first introduced by 
Longuet-Higgins and Steward (1964), is utilized to express the change in momentum due 
to waves. Thornton and Guza (1983) using linear theory and the energy flux equation(2) 


showed that the wave induced radiation stress gradient mav be expressed as 


eS. sin, 2 sin 0, 
Seep eer [14] 
x 2 & 














Neglecting the surface wind stress in the alongshore direction and substituting equation 


(12), the alongshore momentum equation becomes: 


sin 6, ase b 
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The tangential bottom stress t4 for small angle of incidence and weak mean 


longshore current can be approximated ( Longuet-Higgins (1970a,b) ) : 


u 








c= peu lV [16] 


in which c, is the bottom shear stress coefficient, | is the mean water particle velocity 
magnitude and V is the mean longshore current. Thornton and Guza (1983) calculated 
the mean water particle velocity magnitude as : 

mre 
or a7 


u 








= |99 


The fateral turbulent mixing term S’,, is the vertically integrated Reynold’s stress. 
Since the functional form of this stress is not known, it is usually parameterized utihzing 


an eddy viscosity coefficient v : 





’ i, pty. h dV 8 
Syn = see oa pil [18] 
The eddy viscosity formulation by Longuet-Higgins (1970b), which spatially decrease 


from offshore to onshore, is used in the model : 
v= Vixlogh [19] 


where N 1s an adjustable coefficient and based on dimensional analysis has an approxi- 
mate range of O< VN < 0.016 

If the lateral transfer of momentum is not negligible and is to be included, then the 
alongshore momentum balance equation is a second order ordinary differential equation. 


In this case, the gradient of this term is given by: 


aS v oie > AV 
aS pNvg Be al 





dx ax dx 


Taking the derivatives, 


es 3. =. dh av n air 
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Substituting back into the steady state alongshore momentum equation, we have: 
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Grouping the terms. we end up with the following second order differential equation. 
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[23] 


C. APPLICATION OF THE MODEL 
1. Wave Height Model 

To apply the wave height model. the the parameters B and y are combined, re- 
ducing the degree of freedom to one. The B parameter describes the amount of foam on 
the face of a breaking wave i.e. the intensity of wave breaking. y is a proportionality 
constant that relates the rms breaking wave height to the water depth and is slope de- 
pendent. The combined parameter is designed to reduce the number of free parameters 
to tune the model. The term in brackets of equation 12 in the bore dissipation equation 
1S a Correction term to the original bore dissipation formulation. Therefore when the B 
and ; parameters were combined the ; term in the brackets was assigned a typical value 


ae B 
of 0.40, The combined parameter BG is equal to => . 


é 


The model generated wave heights and data were compared by finding the BG 


value that results in the minimum root mean square error 


(ee 1 model ea 
—_ we [26] 
ae - fi Ye E ay 


where /1/,,,4, 18 the model generated wave height, /7,,,. is the measured wave height and 
n is the number of data points. The model was iterated over the combined parameter 
BG and the rms error was computed. The BG parameter was varied from low to high 
values to observe the model behavior and stability. Based on the minimum rms error, 
the best fit was obtained. Examples from each data set for the model fit to the wave 
height data are shown in Figures 9 to 14. 

2. Longshore Current Model 

The longshore current model was compared with the measured longshore cur- 
rents for planar beaches and reasonable results were obtained with and without turbu- 
lent mixing. Similar comparision for barred beaches showed that the model fails to 
predict the longshore current over the shore side of the bar when turbulent mixing is not 
included. In order to improve the results. turbulent mixing parameterized by an eddy 
viscosity was included in the model. which increased the degree of longshore current 
equation to the second. 

This second order differential equation (23) was solved numerically using a 
centered difference scheme. The equation was put in a finite difference form and a linear 
system of equations were obtained for the interior of the spatial domain. At the 
shoreline, a zero current boundary condition was assumed and the very offshore value 
measured in the data was used as the offshore boundary condition. The system of 
equations was put in a tridiagonal matrix. The tridiagonal matrix was solved by a con- 
venient tridiagonal solver by Maron (1987). A mid value of 0.008 was assigned to the 
lateral mixing coefficient, N , which has an approximate range of 0.0 to 0.016 and the 
model was iterated over the bed shear coefficient until the error between the model 
generated current and the data is minimized. Then with the best shear stress coefficient 
obtained the model was iterated over N from 0.0 to .030 until the rms error was reduced. 
Examples of model fits of Superduck and Santa Barbara data are shown in Figure 15 and 


16 respectively. 


Il. DATA AND METHODS 


A. WAVE HEIGHT TRANSFORMATION 

To test and modifv the wave transformation model, 74 data sets were used. The 
data consisted of both laboratory and field data. These data provided many different 
wave conditions from the long period, low steepness Pacific swell to steeper Atlantic and 
North Sea storm waves. The experiments were conducted over various bottom profiles 
including both planar and barred beaches. The data sets are described below. The data 
sets are named after either the site of the data or the person that acquired the data. In 
the tables that describe the data, tanf is the mean heach slope, fis the peak frequency 
in hertz, /,,, is the root mean square wave height, zZ is the deep water wave steepness, 
€, is the deep water surf similaritv parameter, BG is the combined model parameter and 
é,n; 1s the rms error percentage between the model generated wave height and the data. 


The deep water surf similarity parameter 


ee ee [27] 
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was calculated for each data set in which //. is the deep water wave height, L, is the deep 
water wave length and tanf is the mean beach slope (Batyes 1974). The deep water 
wave height was calculated from the measured wave height in shallower water using the 
linear wave theory shoaling coefficient corresponding to the the peak frequency. Simi- 
larly. the deep water wave length was calculated from the peak frequency using linear 
wave theory. 
1. Torrey Pines ( NSTS )Data 

Torrey Pines wave height data were acquired as part of the Nearshore Sediment 
Transport Study (NSTS! at Torrey Pines Beach, CA, during November 1978. The beach 
Was gently sloping. near planar and homogeneous in alongshore direction. Shadowing 
effects of the offshore islands and refraction limits the angle of wave incidence in the surf 
zone to a narrow band of less than + 15 degrees and usually much less. The winds 


during the experiment were light and variable during the experiment. An extensive array 


of instruments was deployed to study the nearshore wave dynamics. The sensors used 
in this study were 11 Marsh-McBirney electromagnetic current meters, four Stathem 
pressure transducers and four resistance wire wave staff. The wave heights were deter- 
mined by the zero-up crossing method, in which the wave height is defined as the dif- 
ference between maximum and minimum within two consecutive zero-up crossings. The 
data were filtered to include only the wind wave band of frequencies between 0.05 and 
0.3 Hz. The data are described in Table 1. Details of the experiment are given in 
Thornton and Guza (1983). 


Table 1. TORREY PINES WAVE HEIGHT DATA 
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2. Santa Barbara (NSTS) Data 
The Santa Barbara data (Table 2) were acquired during the second NSTS ex- 
periment conducted at Leadbetter Beach from 30 January to 23 February 1980. The 
site was selected because the waves are homogeneous in the alongshore direction and the 
bottom contours are relatively straight and parallel. The Leadbetter Beach is composed 
of wellsorted fine to medium size sand. The mean nearshore slope varied between .03 
to .06 during the experiment. The shoreline has an unusual east-west orientation along 
a predominantly north-south coast. The open ocean waves are limited to a narrow 
window of approach -9 degrees to +9 degrees centered at 249 degrees. Therefore the 
waves from the open ocean were directionally narrow banded swell arriving at relatively 
large angles of about 10 degrees to the shoreline. Instrumentation included 24 
Marsh-McBirney electromagnetic current meters and 14 pressure sensors. The exper- 
iment is described in detail in Thornton and Guza (1986). 
3. Batties and Stive (1985) Data 
The first three cases ( Table 3 ) are lab data that were acquired in a wave flume 


using mechanically generated random waves. The bottom profiles include a schematized 


Table 2. SANTA BARBARA WAVE HEIGHT DATA 
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bar profile and a plane slope in concrete (Case 5 and Case 10). Case 15 has a profile 
with bar and trough on sand. All of the lab wave height data were measured by means 
of parallel wire conductivity wave gages. 

Case 18 data is field data acquired on a beach near Egmond City, The 
Netherlands. The field site was chosen to have more or less statistically uniform condi- 
tions alongshore and normal wave incidence. These field data were acquired under 
stormy conditions during which the incident rms wave height was 2.8 meters. It was 
noted that the wind influence was negligible. The //,,, was calculated from the surface 
elevation time series assuming that the surface elevation is a narrow band Gaussian 
process, Which appears to be in good agreement with observations ( Thornton and Guza 
1983). The rms wave height is related *o the standard deviation, ¢, of the surface ele- 


vation by 


Hims = S86 [28] 


Table 3. BATTJES AND STIVE WAVE HEIGHT DATA 
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4. Superduck (1986) Data 

The Superduck data (Table 4) were acquired during the one month long ex- 
periment held by the U.S. Army Corps of Engineers at the Field Research Facility in 
Duck, N.C. in October 1986. The experiment was called Superduck. The site of the 
experiment is along an extensive barrier island formation known as North Carolina’s 
Outer Banks. There are no major structures along the shoreline. The beach has a mean 
foreshore slope of 1:10, usually a single bar formation and a mean slope of 1:100 off- 
shore the bar. Extensive bathymetric data were acquired using the Coastal Research 
Amphibious Buggy (CRAB) so that the wave height and longshore currents could be 
processed with accurate and detailed bathymetry data. The waves and currents were 
measured using a pressure sensor and a Vertical array of three Marsh-McBirney current 
meters mounted on a moveable sled. The details of the data can be found in Whitford 
and Thornton (1989). 


Table 4. SUPERDUCK WAVE HEIGHT DATA 


DATA lem €rms 
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S. Vincent (1984) Data 
The Vincent data (Table 5) were acquired from a series of experiments done at 
the Coastal Engineering Research Center (CERC) in a wave tank of 44 meters long, 0.6 
meters wide and a water depth of 0.6 meters. At one end of the tank was a smooth 
concrete ramp with a slope of 1:30. A random wave generator was used to generate 
waves. The wave height data were measured using nine resistance wave gauges (Vincent, 
1984). 
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Table 5. VINCENT WAVE HEIGHT DATA 
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Figure 3. Setup In CLERC Wave Tank: Ten gages were placed three of them 
were located over the horizontal bottom and seven of them over the 1:30 


slope. 


6. CERC (1983) Data 
Wave height data (Tables 6 and 7) were acquired in Coastal Engineering Re- 
search Center (CERC)’s wave tank. The tank was 0.9 meters high, 0.4 meters wide and 
45.7 meters long. The setup of this tank was designed to obtain data on the shoaling 
and breaking random waves. Ten parallel resistance wire wave gages were used. Three 
of them were located over the horizontal bottom of the tank and seven of them over the 


1:30 concrete slope ( Figure 3 ). 


Table 6. CERC WAVE HEIGHT DATA ( PART I) 
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Table 7, CERC WAVE HEIGHT DATA ( PART II) 
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B. LONGSHORE CURRENTS 

The longshore current model described in Chapter 2 is applied to the Santa Barbara 
and Superduck data sets. These are the only high qualitv longshore current data avail- 
able acquired on a homogeneous alongshore bathymetry with adequate measure of di- 
rectional wave forcing. The Santa Barbara beach profiles were almost planar while the 
Superduck profiles were barred. Of the combined two months of data collection, five 
cases from Santa Barbara and four cases from Superduck were selected for comparision. 

Since the longshore currents are induced by the wind waves, simultaneous meas- 
urements of the wave field and longshore current measurements is essential. This re- 
quirement was satisfied in the data sets. In the Tables 8 and 9 describing the longshore 
currents, c, is the bed shear stress coefficient, N is the lateral mixing coefficient, fis the 


: : die 
frequency in hertz, //,,, is the the rms wave height, —— is the deep water wave steepness, 


rms 


* . . . 6 
€, is the deep water surf similarity parameter, e,,, 15 the rms error percentage between 


rms 


the model generated longshore current and the data. 


Table 8.5 SANTA BARBARA LONGSHORE CURRENT DATA 
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Table 9. SUPERDUCK LONGSHORE CURRENT DATA 
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IV. DISCUSSION 


The wave transformation model was applied to 74 different data sets. The very 
offshore wave height and direction measured were used as an input to the model. The 
model is not verv sensitive to the exact orientation of waves to the shoreline since the 
group speed computations involves the cosine of the angle of wave incidence and the 
cosine of a small angle is not significantly different from one. Therefore exact wave 
Orientation 1s not required and the error due to wave orientation is expected to be small. 

The data sets and the profiles were compiled and the model was iterated over the 
BG parameter to reduce the root mean square error percentages in the least square 
sense. A fully developed bore would be expected to have a B value of | and since the 
typical proportionality constant ; for Rayleigh distribution is about 0.4, the combined 
parameter BG for a fully developped bore would be equal to 6.25. The statistics of the 
BG parameter for the whole 74 data sets were computed. A mean of BG parameter of 
4 with a standard deviation of 2.9 was obtained for all the data sets. 

Model sensitivity to the choice of the parameter was investigated by Gill (1985). 
He varied the y parameter from 0.35 to 0.55 and found that the error increase was less 
than §-10%o. In the same manner, the B parameter was iterated from 0.9 to 1.7 and it 
was found that the increase in rms error was less than 10-12%. Similar type of investi- 
gations were done by Thornton and Guza (1983) for Torrey Pines data set and it was 
found that varving the B parameter + 25°o resulted in an increase of rms error less than 
10° o. 

The combined parameter BG is plotted against the rms error in Figure 5, and it 1s 
seen that changing the BG parameter one standard deviation either side of the mean 
causes an increase in rms error of less than 9% . Based on these investigations it is 
evident that the model is not overly sensitive to the exact choice of parameters. 

Relatively high values of BG parameter are required to minimize the rms error for 
Torrev Pines data. The Santa Barbara data had almost the same deep water steepness 
but didn’t require as high values of BG. The difference is hypothesized due to the defi- 
ciencies of linear bore theory. For spilling type breaker the turbulence is confined to a 
surface layer, basically between the crest and trough of a wave. For plunging or surging 
type breaking the turbulence is not confined between the crest and trough, but pene- 


trates much deeper. The actual depth h. concerning the volume discharge computation 
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Figure 4. Variation of B with Deep Water Surf Similarity Parameter; The B 
parameter shows a linear trend with the deep water surf similarity pa- 
raincter, €,. The dashed hne is the hnear regression line. Notice that the 
B parameter values increase as the breaker type shifts from spilling to 


plunging (¢, > 0.4). 


will be different from the mean water level as expressed by the linear bore theory. It is 
evident that there is deviation from the bore formulation. Although the incident wave 
parameters such as decp water steepness and frequency for Torrey Pines and Santa 
Barbara data sets are about the same, the beach slope of Santa Barbara is about twice 
that of Torrey Pines. This result implies that the wave breaking is a function of both 


Lie ; eee 
deep water Wave steepness (a ) and beach stope (tanf). Therefore the wave dissipation 
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Figure §. Variation In Rms Error Due to Change In BG Parameter: Changing 
BG (BETAGAMMA) parameter one standard deviation on either side 


of the mean docs not increase the error greater than 9%. 


rate 1s lower over mild slopes compared to steep slopes. The combination of plunging 
type breaker and the beach slope effect is thought to be responsible for high values of 
BG for Torrey Pines data set. This idea is confirmed when the BG paraineter obtained 


for all cases plotted versus deep water steepness. As seen in (Tigure 6), as the deep 


2| 


water steepness becomes less than 0.02 the BG parameter tends to increase exponen- 
uallv. For deep water steepness greater than 0.02 the BG parameter stavs invariant at 
a value of about 2.9. The curve fitted to the individual points has the form of 
oo - i where @ is equal to 2.25 and f is equal to 0.0085. 

Using the deep water surf similarity parameter the breaker type is classified as fol- 
lows (Battjes, 1974) : 


&,<0.4 Spilling 
04=25=70 Plunging 
6. 2.0 Surging or Collapsing 


Most of the data sets tabulated in chapter 2 fall in the spilling range with a few of them 
being in the plunging range. There were no surging or collapsing type breaker. Torrey 
Pines and Santa Barbara (NSTS) data differ from the other data due to their low steep- 
ness, about one order of magnitude smaller than the other data sets. Asa result, the 
breaker type based on the surf similarity parameter by Battjes (1974) suggest plunging 
for all cases of the NSTS data sets. Breaking intensity increases from spilling to plunging 
to collapsing breakers as &, increases. Therefore it is expected B and hence BG will be 
a function of €, and increase as €, increases. 

The BG parameters obtained for all data sets are plotted versus deep water surf pa- 
rameter and an exponential trend was observed. A curve in the form of y = ae** is fitted. 
Excluding one outlving data point seen in the right side of Figure 7 , it is found that 
#=1.303 and f= 2.4 fitted the data well. This curve acts in a very sunilar manner to the 
breaking wave height curves in the Shore Protection Manuel, U.S. Army Corps of En- 
gineers (1984) in this exponentially increasing region of the curve. Utilizing this func- 
tion, the model parameter can be predicted. The deep water surf parameter can be 
computed for a particular beach of interest and inputting this value to the fitted curve 
function, the BG parameter can be predicted. 

The BG parameter was plotted versus the incident rms wave height. No correlation 
was evident. Therefore it was concluded that the model parameter is not strongly de- 
pendent on incident wave height by itself Running all the data sets, it was observed that 
the model is capable of handling a small scale wave height of © cm lab data as well as a 


278 cm large wave height measured during a stormy weather. 
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Figure 6. BG Versus Deep Water Steepness: The sohd curve represents the 
curve fitted to BG Mayues versus deep water steepness for all data sets in 


the form of y=o+->. 
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Figure 7. BG Versus Deep Water Surf Parameter: The solid curve represents 
the curve fitted to BG values versus deep water surf parameter for all 


data sets in the form of y = ae". 


WATER DEPTH(CM) 


Figure 8. 


RMS WAVE IIEIGHT AND DEPTIT PROFILE 
TORREY PINES NOV 20 1978 BG=12.2 


a 
oO 
— 
ee 
Be 
2 
re) 
joo} 
f 
> 
< 
= 
n 
= 
oe) 


9800 14700 
Distance Offshore (Cm) 





Rms Wave Height and Depth Profile of Sample Torrey Pines 
Data: Model predicted wave height versus November 20 1978 data. 
The solid line is the mode! predicted wave height while the individual 
points are the measured wave heights. The second plot shows the depth 


profile. 
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Figure 9. Rims Wave Height and Depth Profile of Sample Santa Barbara 
Data: Model predicted wave height versus February 03, 1980 data. 
The solid line is the model predicted wave height while the individual 
points are the measured wave heights. ‘Ihe second plot shows the actual 


depth profile. 


Sallenger ez al. (1985) investigated the dependence of wave height on local slope and 
depth. They found that y linearly increases with bottom slope in the form of: 


y = 3.2 tan f+ 0.30 [29] 


They also concluded that the y coefficient was independent of the deep water wave 
steepness. The value of y was computed for all the data sets using the above relation. 
The coefficients obtained based on this relation ranged from 0.33 to 0.48 with a mean 
of 0.40. Their conclusion that y was independent of the deep water wave steepness was 
found also to be valid for the present data sets tabulated in chapter 3. The B values were 
computed by substituting the calculated y values into the combined BG parameter. The 
computed B parameter values in this manner ranged from 0.46 to 1.24 with a mean of 
0.82 and a standard deviation of 0.17. This result shows that B is assigned reasonable 
values and agrees with the physics of the problem. These B parameter values showed 
an increasing trend with the deep water surf similarity parameter (see Figure 4). As the 
breaker type shifts from spilling to plunging (€, > 0.4), the B parameter values increase. 

The longshore current model was run first neglecting turbulent mixing. Relatively 
high rms error percentages in the least square sense were obtained (See Figure 15 ). 
No matter how much c, was adjusted, it was not possible to obtain longshore current 
over the shore side of the bar where the waves reform and are no more breaking. Based 
on this result, it was decided to include turbulent mixing in the model. McDougal e 
al. (1986) investigated the influence of lateral mixing on longshore currents. They tested 
a variety of latera! mixing models and concluded that the longshore current profile is 
insensitive to the form of the mixing model. Therefore it was decided to use the lateral 
mixing model of Longuet-Higgins (1970b) which has the form of vy = Ax/gh . Inclusion 
of the turbulent mixing increased the degree of the longshore current equation to the 
second. More reasonable results are obtained by including the turbulent mixing term 
(Figure 15). 
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Figure 10. 
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Rins Wave Height and Depth Profile of Sample Battjes and Stive 
Data: Model predicted wave height versus Case 18 data. The solid 
line is the model predicted wave height while the individual points are 
the measured wave heights. The second plot shows the actual depth 


profile. 
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Figure 11. Rms Wave Height and Depth Profile of Sample Superduck 


Data: Model predicted wave height versus October 16, 1986 data. 
The solid line is the model predicted wave height while the individual 
points are the measured wave heights. The second plot shows the 


depth profile. 
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Figure 12. Rims Wave Height and Depth Profile of Sample Vincent Data: Model 
predicted wave height versus Case 1509 data. The solid line is the 
model predicted wave height while the individual points are the meas- 
ured wave heights. The second plot shows the constant slope bottom 


profile. 
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Figure 13. Rms Wave Height and Depth Profile of Sample CERC Data Part 
(I): Model predicted wave height versus Case 44 data. The solid line 
is the model predicted wave height while the individual points are the 


measured wave heights. The second plot shows the depth profile. 
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Figure 14. Rms Wave Height and Depth Profile of Sample CERC Data Part 
(Il): Model predicted wave height versus Case 68 data. The solid line 
is the model predicted wave height while the individual points are the 


measured wave heights. The second plot shows the depth profile. 
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Figure 15. Longshore Current and Depth Profile of Sample Superduck Data: The 
dotted line indicates the model generated longshore current without 
turbulent mixing. The solid Jine ts the model predicted current with 
turbulent mixing and the individual points are the measured wave 
heights. The second plot shows the depth profile. Notice that the ve- 
locitics on shore side of the bar significantly change depending on 


inclusion or exemption of the turbulent mixing. 
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Figure 16. Longshore Current and Depth Profile of Sample Santa Barbara 
Data: ‘the dotted line indicates the model generated longshore current 
without turbulent mixing. The solid line is the model predicted current 
with turbulent mixing and the individual pomts are the measured wave 


heights. The second plot shows the depth profile. 


In the Santa Barbara data, the model generates a maximum longshore current closer 
to the shoreline compared with the maximum current measured (see Figure 16). This 
result was regardless of the inclusion or exemption of the turbulent mixing. Therefore, 


the inclusion of turbulent mixing to the model did not cause significiant improvement 


for Santa Barbara longshore current data. The model comparision was found insensitive 
to the N value. Based on this result, turbulent mixing may be excluded in the model for 
planar beaches. As disscussed earlier this is not true for barred beaches. This would 
reduce the number of model parameters to be specified and also since the longshore 
current equation will stay in the first order, some cpu time will be saved. Planar beaches 
are more appropriate for Amphibious Warfare compared to barred beaches. Landings 
will not be made on barred beaches unless necessary. If a barred beach is to be chosen 
for landing. then the turbulent mixing should be included. 

The mean value of the bed shear stress coefficient, c, for the present longshore cur- 
rent data is 0.006 and the mean lateral mixing coefficient is 0.0055. Since there is not 
enough longshore current data to obtain predictive relations for the model parameters, 
these mean values are suggested for use in model implementation. Longshore current 
magnitude is directly proportional to c, It should be recognized that, since the optimal 
c, Values varied between 0.0025-0.009, the expected current magnitude can have 
significiant error using the suggested values. For generality, since it does not matter to 
include or exclude the turbulent mixing term for planar beaches, the model can run with 
these c, and N values regardless of the beach profile. These values should be checked 


as more longshore current data becomes available. 


V. SUMMARY AND CONCLUSIONS 


In this study, the wave height transformation and longshore current models of 
Thornton and Guza (1983, 1986) are applied to a varietv of wave conditions over planar 
and barred beaches for the purposes of testing and evaluation. 74 data sets of wave 
height transformation and 9 data sets of longshore currents are utilized. The model re- 
sults are optimized in the least square sense to obtain the best fit to the data. The model 
coefficients in the wave transformation model, y a measure of saturated wave breaking 
and B a measure of breaking wave intensity, are combined into a single BG parameter. 
Based on data analysis, the combined BG parameter in the wave transformation model 
is found to be correlated with the deep water surf parameter. Therefore, this correlation 
can be used to predict the single model BG parameter in future model applications. 

Testing the model using field and laboratory data shows that the wave transforma- 
tion model is highly robust in describing the wave height distribution inside the surf 
zone. Varying the combined BG parameter one standard deviation from the mean of the 
74 data sets, did not increase the rms error greater than 9%. Therefore it is concluded 
that model 1s not overly sensitive to the parameters. 

The » parameter was calculated using the relationship by Sallenger and Holman 
(1985) and the B values were calculated from the BG parameter. The mean propor- 
tionality coefficient y for the whole 74 data sets was 0.4. The B parameter, which de- 
scribes the intensity of wave breaking. appears also to be correlated with the deep water 
surf parameter. This is expected since the deep water surf parameter is used to predict 
breaker type. and B is a measure of breaking wave intensity. 

lor longshore current computations the turbulent mixing needs to be included over 
barred beaches, but for planar beaches its inclusion has an insignificant effect on the 
results. Therefore it is suggested the turbulent mixing be included in a general model for 
consistency. At present there is not a means for predicting the appropriate model pa- 
rameters for longshore current computations. Therefore it is suggested the longshore 
current model be run with the mean values obtained for the bed shear stress and lateral 
mixing coefficients, 0.006 and 0.0055 respectively, regardless of the beach profile. 

The model is operationally handy since it is only initialized by the offshore wave 


height, direction and period, which can be observed aboard the ship. The construction 


of the model is quite simple, and does not require much computation time. Therefore 


the model can run on a personal computer with almost instant response. 
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